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PRESSURE GRADIENT EFFECTS ON HEAT TRANSFER TO REUSABLE 
SURFACE INSULATION TILE-ARRAY GAPS* 

David A, Throckmorton 
Langley Research Center 

SUMMARY 

An experimental Investigation was performed to determine the effect of pressure 
gradient on the heat transfer within space shuttle reusable surface insulation (RSI) tile- 
array gaps under thick, turbulent boundary-layer conditions. Heat-transfer and pres- 
sure measurements were obtained on a curved array of full-scale simulated RSI tiles in 
a tunnel-wall boundary layer at a nominal free-stream Mach number of 10.3 and free- 
stream Reynolds numbers of 1.6, 3.3, and 6.1 x 10® per meter. Transverse pressure 
gradients of varying degree were induced over the model surface by rotating the curved 
array with respect to the flow. Definition of the tunnel-wall boundary -layer flow was 
obtained by measurement of boundary -layer pitot pressure profiles, wall pressure, and 
heat transfer. 

Flat -plate heat-transfer data were correlated and a method was derived for pre- 
diction of heat transfer to a smooth curved surface (i.e., no gaps) in the highly three- 
dimensional tunnel -wall boundary -layer flow. Pressure on the floor of the RSI tile- 
array gap followed the trends of the external surface pressure. Heat transfer to the 
surface immediately downstream of a transverse gap is higher than that for a smooth 
surface at the same location. Heating to the wall of a transverse gap, and immediately 
downstream of it, at its intersection with a longitudinal gap is significantly greater than 
that for the simple transverse gap. No systematic effect of transverse pressure gra- 
dient on heating was observed. Simulation of full-scale space shuttle vehicle boundary- 
layer thickness and pressure -gradient level was good. 

INTRODUCTION 

The space shuttle orbiter is being designed for an operational life in excess of 
100 reentry missions with minimal refurbishment between those missions. In order to 

Part of the information presented herein was included in a thesis entitled "An 
Experimental Investigation of the Heat Transfer to Reusable Surface Insulation Tile- 
Array Gaps in a Turbulent Boundary Layer With Pressure Gradient," submitted in par- 
tial fulfillment of the requirements for the degree of Master of Science, The George 
Washit^ton University, February 1975. 



meet this requirement, the vehicle thermal protection system (TPS) will be a surface 
covering of a nonmetallic, low-density refractory oxide. This material, referred to as 
reusable surface insulation (RSI), must be capable of withstanding, without degradation, 
repeated exposure to the harsh reentry environment, while insulating the structure of the 
vehicle from surface temperatures in excess of 1500 K. The material will be attached to 
the vehicle surface in a "brick-like" array of square tiles (15.25 by 15.25 cm) which vary 
in thickness from approximately 1 to 10 centimeters according to the Intensity of the local 
heating. Small gaps between tiles will allow for thermal expansion of the tile material. 
Interference heat transfer to the tile exterior and gap wall surfaces is of major concern 
to the TPS designer. The presence of the gaps may result in increased boundary -layer 
turbulence and attendant increased surface heating. Heating levels within the gaps, which 
would be expected to be substantially lower than surface values, may approach surface 
values due to flow reattachment phenomena. In addition, radiation blockage within the 
gaps may produce extreme temperatures even at low heating levels, and the shortened 
heat paths (i.e., from gap wall instead of tile exterior surface) may result in excessive 
bond-line temperatures. 

Effective design of the TPS requires a sound knowledge of the aerodynamic heating 
environment to which the RSI tiles will be subjected. This knowledge must include an 
accurate definition of the heat-transfer distribution within the tile gaps and a good under- 
standing of how this distribution is affected by boundary -layer state (laminar/turbulent), 
boundary -layer thickness, flow angularity, gap width, tile edge radius, tile stacking 
arrangement, and other parameters. 

As a part of the space shuttle development program, an experimental effort has 
been focused on the shuttle -related gap heating problems. Evaluation of TPS tile material 
thermal performance (ref. 1) and definition of the gap heating environment are important 
elements in this effort. Johnson (ref. 2) studied the effects of gap width and boundary- 
layer thickness on TPS gap heating for turbulent flow over a simulated tile array at 
Mach 8. This work revealed a potential heating problem area at the intersection of 
streamwise and transverse gaps. The present author obtained data on a simulated tile 
array in a turbulent boundary layer to investigate the effects of gap width, flow angularity, 
and tile surface mismatch on tile gap heating. (See ref. 3.) Foster et al. (ref. 4) investi- 
gated the effect of similar parameters for single gaps and gap intersections in laminar 
flow. A compilation and analysis of these and other recent gap heating data may be found 
in reference 5. 

Because of entry angles of attack and the complex curvatures of the orbiter external 
surface, boundary -layer flows over much of the flight vehicle will be strongly influenced 
by surface pressure gradients. Each of the recent studies of TPS gap heating phenomena, 
however, has been conducted on sharp-leading-edge flat plates or in tunnel-wall boundary 
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layers under conditions of zero pressure gradient. The present investigation was under- 
taken to assess possible effects of pressure gradient on the heat transfer to RSI tile -array 
gaps in a thick, turbulent boundary layer. 

Heat-transfer and pressure-measurement tests were made on a curved array of 
full-scale simulated RSI tiles submerged in a thick, turbulent, tunnel -sidewall boundary 
layer. Transverse pressure gradients of varying degree were induced over the model 
surface by rotating the curved array with respect to the flow direction. This enabled 
the study of RSI tile-array gap heating as affected by pressure gradient. The tile -array 
model was tested both with gaps present, and with the gaps filled and smoothed to pro- 
vide smooth-surface reference data. Heat-transfer and pressure measurements were 
also made on a flat plate mounted in the tunnel sidewall. These data, along with meas- 
ured pitot pressure profiles, provided a definition of the characteristics of the three- 
dimensional boundary -layer flow in which the tile-array tests were conducted. 

SYMBOLS 

A constant in equation (7), cm"^ 

B constant in equation (7) 

Cp constant pressure specific heat of air, j/kg-K 

Cp,m specific heat of model material, J/kg-K 

d depth into gap, cm (see fig. 5) 

h heat-transfer coefficient, W/m^-K 

k constant in equation (14) 

L full-scale tile dimension, cm 

M Mach number 

Ng^ Stanton number, h/^^UgCp 

P pressure, N/m^ 

^t,2 total pressure behind a normal shock in the free stream, N/m^ 
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q 

R 




r 

S 

T 

t 


u 

X,Y,Z 

x,y,z 

^max 

^surf 

a 


7 

f 

5eff 

d 


heat-transfer rate, W/m^ 
universal gas constant, mVsec^-K 

Reynolds number based on wall conditions and boundary -layer momentum 
thickness, p^Ug0/fXe 

free -stream Reynolds number, p^u^y/p^, m"^ 
recovery factor 

surface dimension defined in figure 30, m 
temperature, K 
time, sec 
velocity, m/sec 

coordinate system defined in figure 10 

coordinate system defined in figures 3 and 10 

maximum vehicle planform dimension defined in figure 30, m 

z -dimension of local curved model surface, cm 

array rotation angle, deg (see f^. 10) 

ratio of specific heats of air 

boundary -layer displacement thickness, cm 

"effective" boundary-layer displacement thickness defined in equation (14), cm 
boundary -layer momentum thickness, cm 
model material thickness, m 
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H viscosity, N-sec/m2 

p density of air, kg/m^ 

Pj^ model material density, kg/m^ 

0 angle between surface tangent plane and flow velocity vector, deg 

Subscripts: 

aw adiabatic wall 

e boundary -layer edge 

fp flat plate 

sm surface of the smooth, curved array 

t total 

w wall 

°o free stream 


APPARATUS AND TESTS 
Facility 

The experimental results presented herein were obtained in the Langley continuous 
flow hypersonic tunnel. This facility, which has a 78.74-cm square test section, operates 
at a nominal free-stream Mach number of 10.3 over a free-stream Reynolds number range 
of 1.5 to 8.2 X 10® per meter using air as the test gas, and may be operated in either a 
blowdown or continuous, closed-circuit mode. To prevent liquefaction, the air is heated 
by use of electrical resistance tube bundles. The tunnel-throat, expansion, and diffuser 
sections are all water cooled. A photograph of the facility test section is presented in 
figure 1. 

For these tests, the models were mounted on the model injection mechanism which 
is shown adjacent to the test section in figure 2. This device allows a model to be iso- 
lated from the hypersonic air stream for model cooling or configuration changes while 
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the tunnel is operating. The mechanism also provides rapid injection of a model into 
the hypersonic airstream. 


Models 

Pressure -gradient model .- The RSI tile pressure-gradient model was a curved 
array of simulated full-scale RSI tiles fabricated of Inconel sheet with a nominal thick- 
ness of 0.0483 cm. Surface curvature was generated from a r^ht circular cylinder of 
102.28-cm radius, cut by a plane at a 5° angle to the axis of the cylinder. A schematic 
drawing of the model is shown in f^ure 3. Individual tile size was 14.92 cm square. 

The model was fabricated such that each transverse row of simulated tiles was a 
continuous sheet of material with the longitudinal gaps formed by bending. The transverse 
gaps were fabricated separately and electron-beam welded between tile row sections to 
form the complete tile array. Stress relief in the material following the welding process 
resulted in a "pinching down" of the transverse gaps from the desired width. Resulting 
gap width between tiles was nominally 0.30 cm for longitudinal gaps and 0.20 cm for trans- 
verse gaps. Gap depth was 2.86 cm. 

Spacer plates were also fabricated to maintain model edges flush with the tunnel 
sidewall. Three sets of spacer plates were fabricated to allow testing at model flow 
angularities of 0°, ±22.5°, and ±45°. The model and spacer plates mounted on the facil- 
ity injection system are shown in figure 4. 

The model was instrumented with 67 chromel-alumel thermocouples spotwelded to 
the model back surface at tile-surface and gap locations indicated in figure 5. The model 
was also fitted with 12 static pressure orifices (7 located on the surface of a tile and 
5 located on the floor of a transverse gap) as indicated in figure 5. Precise thermocouple 
and pressure orifice locations are listed in ta.ble I. 

For smooth -surface testing of the model, the gaps were eliminated by filling them 
with twine to a point approximately one gap width below the tile exterior surface. The 
remainder of the gap was filled with plaster of parts which was allowed to dry and then 
smoothed by sanding to match the surrounding surface. This created a smooth model 
with correct curvature for measurement of surface reference data. 

Flat-plate model. - A smooth flat-plate model, used to measure undisturbed wall- 
heat-transfer coefficients, was fabricated from 321 stainless steel sheet with a nominal 
thickness of 0.127 cm. The 50.8-cm square panel fit flush with the injection plate fab- 
ricated for testing of this model. The model and injection plate are shown mounted on 
the injection strut, ready for testing, in figure 6, The flat-plate model was instrumented 
with 11 chromel-alumel thermocouples spotwelded to the model back surface, located 
as shown in figure 7. 
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Boundary-layer probe .- An 11 -tube pitot-probe rake was placed downstream of the 
flat-plate model at four spanwise locations. Individual pitot probes were located at sta- 
tions normal to the wall as shown in figure 8. Tube diameters varied with position rela- 
tive to the wall. The rake is shown installed with the flat-plate model for testing in 
figure 6. 


Instrumentation 

Temperature data were obtained using the chromel-alumel thermocouples with a 
reference junction of 324.8 K. The reference-junction temperature was thermostatically 
controlled to within ±0.55 K. 

Pressures were measured using baratrons mounted on the injection strut imme- 
diately behind the model. The baratron is a capacitance -type transducer, operated in 
conjunction with a signal conditioner to allow measurement of pressure over seven ranges 
from 0 to 68.9 N/m2 to 0 to 68.9 kN/m^. The signal conditioner provides automatic rang- 
ing such that the measurements were obtained on the lowest possible range. Measure- 
ment accuracy is estimated to be within ±7 percent of reading. 

Test Procedures and Conditions 

The transient-calorimeter technique was used to measure heat-transfer rates to 
the surfaces of the thin-skin models. The tests were conducted with the models initially 
at room temperature, isolated from the hypersonic airstream within the injection chamber, 
at a pressure equal to the test-section static pressure. With the hypersonic flow estab- 
lished in the test section, the model was rapidly injected to the test position, flush with 
the tunnel sidewall, and temperature and pressure data were automatically recorded on 
magnetic tape by an analog -to -digital converter at a rate of 20 samples per second. The 
model was exposed to the airstream for an interval of time sufficient to allow pressure 
transducer outputs to "settle out," and then retracted. 

Both the tile-array and flat-plate reference models were tested at total pressures 
within the settling chamber of 2.41, 5.17, and 9.65 MN/m^, corresponding to free-stream 
Reynolds numbers of 1.6, 3.3, and 6.1 x 10® per meter. Additional flat-plate data were 
obtained at a total pressure of 12.07 MN/m^ (Reynolds number of 7.4 x 10® per meter). 

The tile-array model was tested at flow angles a of 0°, ±22.5®, and ±45°. Testing 
at positive and negative flow angles allowed measurement of both temperature and pres- 
sure data over the entire model surface while instrumenting only one side of the plane of 
symmetry for either temperature or pressure as shown in figure 5. By "mirror Imagli^" 
the data obtained at negative flow angles, data for the entire model surface were obtained 
for flow angles of a = 0°, 22,5°, and 45°, 


7 



DATA REDUCTION 


Boundary -Layer Profiles 

Boundary -layer velocity and density profiles were derived from wall boundary- 
layer pitot pressure measurements. These calculations required assumptions of the 
distribution of static pressure and temperature through the boundary layer. As shown 
in figure 9, the measured wall static pressure was significantly higher than the free- 
stream static pressure computed from one -dimensional flow considerations. However, 
pressure measurements in the free stream of the facility (ref. 6) show good agreement 
with the computed free-stream static pressure. Static pressure within the boundary 
layer was assumed to vary linearly between the measured wall and computed free-stream 
values. 

Total temperature within the boundary layer was assumed to vary as the square 
of velocity, that is, 


Tt - Tw ^ /u_\2 (1) 

’^t,e ■ ^w 

which is characteristic of turbulent tunnel-wall boundary layers (ref. 7). Unpublished 
total temperature measurements made by D. H. Crawford in the wall boundary layer of 
this facility show good agreement with the quadratic temperature -velocity relation. 

With the measured pitot pressure and assumed static pressure at each point, the 
Rayleigh pitot equation was applied to calculate local Mach number. Local velocity and 
density were then computed using the assumed temperature distribution and the perfect- 
gas equation of state. 


Heat-Transfer Data 

The test procedure of rapid injection of the isothermal model to the test position 
provided a step input in heat transfer to the thin-skin model. Heat-transfer rates were 
determined by the transient-calorimeter technique of measuring the time rate of change 
of the model skin temperature. For data reduction purposes, the one-half second inter- 
val of temperature data Immediately following model injection was disregarded to allow 
flow conditions to stabilize in the gap. This time is in excess of the required time as 
reported by Nicoll (ref. 8). A quadratic least- squares curve was fit to the subsequent 
4-sec interval of data for each thermocouple. Rates of change of temperature with 
time 9Tw/9t were evaluated analytically from the curve-fit expressions at the initial 
point of each curve fit. Heat-transfer rates were then computed from the expression 
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9Tw 


( 2 ) 


The long interval (4 sec) of temperature -time data used for heat-transfer rate calculation 
allowed measurement of the low heating rates found within the gaps which were not dis- 
cernible when a shorter interval (1 sec) of data was considered. An assessment of con- 
duction effects which result from the long data interval, for representative data from this 
test, indicates a maximum error in computed heating rate of less than 10 percent for the 
"worst case" thermocouple. 

Heat-transfer data are expressed in the form of the heat-transfer coefficient h 
defined as 


h = 


q 


“^aw 


"w 


Adiabatic wall temperature was computed from the relation 


(3) 


r = 



(4) 


where recovery factor r was assumed equal to 0.89 for the turbulent test conditions. 
This relation was applied to both flat-plate and curved-array data. 

RESULTS AND DISCUSSION 
Coordinate Definition 

In the discussion of experimental results which follows, the data are referenced to 
coordinate systems defined in figure 10. All flat-plate and wall boundary -layer data are 
referenced to a stream-oriented coordinate system (X,Y,Z) which is fixed within the flow. 
All curved-panel data, both smooth and with gaps present, are referenced to a coordinate 
system (x,y,z) fixed within the rotating curved array. The origins of both coordinate 
systems are located on the tunnel sidewall at the center of rotation of the tile -array 
model. The coordinate pairs of any point (x,y and X,Y) are related by the rotation 
transformation: 


X = X cos a - y sin a 
Y = X sin ot + y cos a 
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Boundary -Layer Surveys 


Boundary -layer velocity profiles measured on the sidewall center line are pre- 
sented in figure 11 for the four free-stream unit Reynolds numbers at which flat-plate 
data were obtained. The profile shapes are characteristic of a fully developed turbulent 
boundary layer and the profiles tend to "fill out" (boundary -layer thins) with increasing 
unit Reynolds number. Profiles obtained at several transverse locations at a single 
unit Reynolds number are compared in figure 12. A significant change in the boundary - 
layer profile shape as a function of transverse position is observed. This transverse 
variation in boundary -layer profile is an indication of the three-dimensional character 
of the boundary -layer flow in a nozzle of square cross section. The three dimension- 
ality results from the corner interaction of the tunnel sidewall, floor, and ceiling 
boundary layers. 

The measured profile data were integrated to obtain values of the bulk quantities 
of boundary -layer displacement 6* and momentum 9 thicknesses, defined by the 
relations 
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+ 



jPU_\ 

PeV 


dZ 


(5) 



The variations of these quantities with free-stream unit Reynolds number and transverse 
position are presented in figure 13. Momentum thickness decreases, as expected, with 
increasing unit Reynolds number, and shows little transverse variation. Boundary -layer 
displacement thickness does not change s^nificantly with unit Reynolds number, but does, 
however, decrease rapidly with transverse position reflecting the changes in boundary- 
layer profile shape noted in figure 12. 

Flat- Plate Results 

Heat-transfer data .- Measured flat-plate heat-transfer coefficients are shown in 
figure 14 for the full range of test conditions. At each Reynolds number, the heat trans- 
fer increases in the transverse direction, primarily as a result of the thinning of the 
boimdary layer, and is symmetric about the center line. The streamwise distribution of 
heat transfer shows a small heating decrease characteristic of a fully developed, thick- 
ening boundary layer. 
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Pressure data .- The flat-plate pressure data (fig. 15) show a transverse distribu- 
tion similar to that of the heat-transfer data. Streamwise, however, the pressure 
decreases to a minimum and increases in the downstream (+X) direction. This behavior 
is attributed to the fact that the test-section area of this facility does not represent a com- 
plete expansion of the flow from the tunnel throat. The tunnel nozzle expands the flow to 
a point approximately 1 meter upstream of the test-section center, at which the flat, par- 
allel walls of the test section begin. The intersection of the expansion- section and test- 
section walls then constitutes a compression corner for the wall boundary-layer flow. 

The data presented herein were obtained in the compression region downstream of this 
corner. 

Flat -plate data correlation . - In order to understand the mechanisms controlling the 
wall boundary-layer heat transfer, the heating data of figure 14 were correlated with 
parameters which could be expected to influence heat transfer. Figure 16 presents flat- 
plate heating data expressed in the form of Stanton number Ngj plotted as a function of 
Reynolds number based on wall conditions and boundary -layer momentum thickness q . 
The heat-transfer data points were Inferred from the data of figure 14 for the transverse 
locations at which boundary -layer profiles were obtained. 

At each transverse location, Stanton number and Reynolds number correlate the data 
over the range of free -stream flow conditions with Stanton number varying as the wall 
Reynolds number R^ g to the -0.07 power as indicated in figure 16. Failure to corre- 
late the data for the three transverse locations is not surprising, as the Stanton-Reynolds 
number correlation relates changes in streamwise variables and has no application to 
transverse flow phenomena. The upstream histories of the boundary layers affecting 
each transverse location are unique and no simple boundary -layer parameter relates the 
transverse influence of boundary -layer characteristics. The single parameter which 
varies significantly in the transverse plane is boundary-layer displacement thickness, 
and, for a given boundary -layer flow, heat transfer is known to decrease with increasing 
boundary -layer thickness. As previously noted, for this tunnel-wall flow, displacement 
thickness does not vary significantly with unit Reynolds number. The effect of transverse- 
boundary -layer -thickness variation on tunnel -sidewall heat transfer is indicated in fig- 
ure 17 where the data of f^ure 16 are presented as a function of displacement thickness. 
The flat-plate heat-transfer data correlate with pertinent flow variables according to the 
relation 


- A6* + B 

where 

A = -0.000277 cm-1 
B = 0.00466 


( 7 ) 
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This expression quantifies the relationship between measured boundary -layer parameters 
and measured wall heat transfer for this tunnel-wall flow. The correlation has no appli- 
cation to any flow other than this tunnel-wall boundary layer; however, it indicates a good 
understanding of this flow and lends confidence to the curved-array data which follow. 


RSI Tile Array — Smooth-Model Data 

Because of the complexity of the wall boundary -layer flow indicated by the flat-plate 
and boundary -layer -probe results, heat-transfer and pressure data were obtained on the 
RSI tile-array model with no gaps present. These data were used to relate measured 
surface and gap heating to an undisturbed surface reference. 

Pressure data.- Smooth-model pressure data and fourth-order least-squares curve 
fits of the data are presented in figure 18 for all test conditions. The data from the ori- 
fice located at y = ±5.08 cm were not considered for computation of the curve-fit 
expressions as data from this transducer were consistently higher than those from the 
other transducers. At an array rotation angle a of 0®, the pressures show the same 
transverse trends demonstrated by the flat-plate data. Flow deflection angle <p (i.c.> the 
angle between the free-stream velocity vector and the plane of tangency of the surface at 
a point) has negligible point-to-point variation between orifices when a = 0°. With 
increases in rotation angle to 22.5° and 45°, flow deflection angles on the upstream por- 
tion of the model increase, those on the downstream portion decrease, and static pres- 
sures vary accordingly.^ 

Heat-transfer data .- Measured smooth-surface reference heat-transfer data are 
presented in figure 19 for the curved-array model. For the zero-rotation case, the flat- 
plate heatii^ data are also presented for comparison. Because of the 5® surface inclina- 
tion of the curved panel, heat transfer to the panel was greater (as shown) than to the flat 
plate. Based on the flat-plate results, the increase in heat transfer with transverse dis- 
tance from the center line was expected. However, at or = 0° the rate of increase in the 
transverse direction was expected to be nearly equal to that for the flat plate since flow 
deflection angle remains constant with transverse position. As figure 19 {a = 0°) shows, 
the transverse rate of Increase of heating to the curved array is significantly greater 
than that for the flat plate. This more rapid increase is hypothesized to be the result of 
an effective "thinning" of the boundary layer due to the protrusion of the model into the 
boundary -layer flow. With a boundary layer which is much thicker than the protruding 
height of the tile -array model, it is not thought that the model will significantly affect 

^Note that the y-coordinate is fixed in the rotating array, while the Y-coordinate 
is fixed in the tunnel sidewall. While y-values are constant for each orifice, the orHice 
locations within the complex boundary -layer flow vary with array rotation angle. There- 
fore the pressure and heat-transfer variations shown in the figures as functions of the 
y-coordinate are a superposition of effects of changing flow deflection a^gle and changes 
in the Y-coordinate position (boundary -layer conditions affecting the point). 
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the outer portions of the boundary layer as would be the case if the characteristic dimen- 
sion of the model was of the same order as the boundary -layer thickness. Rather, the 
boundary -layer edge location remains essentially unchanged from the flat-plate case and 
the boundary -layer thickness is decreased by the protrusion of the curved model into the 
flow. 

Smooth -surface heating data for the a = 22.5° and 45° cases show increased sur- 
face heatii^ with increasing flow deflection angle (y increasing) as expected. The data 
also show, however, increases in heatir^ where the flow deflection angle is decreasing 
(y decreasing). This anomaly is explained by the fact that these heating increases, in 
regions of decreasing flow deflection, are occurring at wall positions where boundary- 
layer thickness is decreasing. The opposing effects of decreasing boundary -layer thick- 
ness and flow deflection angle are dominated by the boundary -layer thinning effect, and 
heating Increases. In the following section, a method will be developed to predict smooth - 
model-surface heat transfer. 

Prediction of smooth-surface heating characteristics .- Consider the correlation of 
flat-plate heating data presented in figure 17 

NstHw" = A6* + B 

and assume that a correlation of this form is valid for the smooth-curved-panel data. 
Therefore, for curved -panel data, 


TJ d0*07 
%^w,0 


^ 0 . 07 \ 
St"w,0/ 


I N Of R, 


fp 


A5eff+ B 
(A6* + B)fp 


( 8 ) 


where is an effective boundary -layer thickness as hypothesized in the previous 

section. By definition. 






R 


w,0 


_^w^ 

Mw 


( 9 ) 


13 



Substituting these expressions in equation (8) and assuming u^g * ^oo,fp’ 


1 

p- 

1 

1 

1 

0.07 

A6*ff+B 

C'Mp 

(^w7'^w)fp 


■ (A6* + B)fp 


( 10 ) 


Assuming a perfect gas, 



and assuming wall-temperature variations are negligible, 


and 


T^-T, 


w,fp 


^w ~ ^w,fp 


Then, 


( 11 ) 


h/P 

PS 

(h/P)fp_ 



A6*ff+ B 
(A6*+ B)jp 


Neglecting the weak dependence on 9 
h ^ A5*ff + B 

hfp \Pfpj (A6*+B)fp 


( 12 ) 


(13) 


Now, estimate 

6*ff = 6* - kzgurf 

where Zg^j-f is the local protrusion of the smooth model surface into the boundary -layer 
flow, and k is an unknown constant. 

If the value of the constant k and the local pressure on the curved array are known, 
the distribution of heat transfer to the surface of the smooth curved panel may be predicted 
using equations (13) and (14) with the measured flat-plate pressure, heat transfer, and 
boundary -layer data. Surface pressure was measured, however, on one transverse ray of 
the smooth curved panel. Figure 20 presents measured heat-transfer data for that ray 
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for a Reynolds number of 3.3 x 10® per meter, along with the predicted heating from equa- 
tions (13) and (14), using the measured pressure data. The value of the constant k which 
provided a "best fit" to the experimental data was 0.5. The predicted heat transfer shows 
excellent agreement with the experimental data. Similar results were observed in the 
data obtained at Reynolds numbers of 1.6 and 6.1 x 10® per meter. 

RSI Tile-Array Results — Gaps Present 

Pressure data .- Measured pressures for the simulated RSI tile array are presented 
in figure 21 for all test conditions. Tile -surface data are shown as open symbols and gap- 
floor data, as solid symbols. Comparison of these data with those of figure 18 indicates 
that the surface pressure distributions with gaps present are essentially identical to those 
obtained with a smooth model. Pressure level within the gaps closely follows that of the 
tile exterior surface. 

Gap heat-transfer data. - Measured distributions of heat transfer^ along the instru- 
mented transverse gap are shown in figures 22 to 24 for each array rotation angle and 
Reynolds number. For the zero-rotation case (fig. 22), at transverse positions off the 
center line (y * 0), no intense reattachment heating occurs at the gap corner, and no sig- 
nificant increase in sirnface heating occurs downstream of the transverse gap, indicating 
negl^ible flow separation over the transverse gap. Note also that off the center line at 
a depth into the gap of less than four gap widths (Ds-symbol), the local heating level is 
less than 2 percent of the undisturbed surface value. The solid symbols in figures 22 
throu^ 29 indicate extremely low heat -transfer -rate data which are of questionable 
accuracy. These data are included for completeness. 

For the zero-rotation case (fig. 22) on the array center line (y = 0) the gap geometry 
consists of a longitudinal gap intersecting the transverse gap. Heat transfer to the trans- 
fer gap at the intersection point and to the tile surface immediately downstream of the gap 
is significantly higher than for a simple transverse gap as exists away from the center 
line. This heating increase at the gap intersection is attributed to a large region of flow 
separation within the longitudinal gap and a strong shear -layer reattachment on the 
forward-facli^ wall of the transverse gap. Heating at the intersection point is depend- 
ent upon Reynolds number. Dunavant and Throckmorton (ref. 9), using gap-intersection 
data from several facilities, have shown that such data may be correlated as a function 
of boundary -layer displacement thickness, streamwise-gap running length, gap width, 
and depth location within the gap. 


^All heat-transfer data discussed in this section are nondimensionalized by the 
measured heat transfer to the smooth- surface model. Surface data are normalized by 
the smooth-model measurement of the same thermocouple; gap-wall data are normalized 
by the smooth-model-surface measurement at the thermocouple location nearest the gap. 
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Rotation of the tile array (figs, 23 and 24) produces increased surface reattachment 
heating immediately downstream of the transverse gap. As in the zero-rotation case, 
the surface reattachment heatii^ is most severe downstream of the longitudinal -gap— 
transverse-gap intersection. Within the gap, array rotation does not significantly affect 
the level of heat transfer, but it does result in a transverse shift of the heating peak on 
the gap wall downstream of the longitudinal gap, toward the low-pressure side of the gap 
intersection. Again, the heating level in this region of flow reattachment is dependent 
upon Reynolds number. 

Figures 25 and 26 present heat-transfer data for the longitudinal gap at each 
Reynolds number as a function of array rotation angle. The significant variations in 
gap heating observed in these plots result from the chaise in orientation of the wall, 
from one which is forward facing to one which is rearward facing, as array rotation 
angle changes. With the exception of certain low Reynolds number data, when the wall 
is rearward facing (a > 0°), the flow appears to separate upstream of the tile gap cor- 
ner, resulting in substantially lower heatii^ to the corner (A-symbol) than to the undis- 
turbed surface. When the wall is forward facing (a < 0°), the flow reattaches in the 
vicinity of the corner resulting in heatii^ equal to or greater than the undisturbed sur- 
face value. 

At a = -45° the longitudinal gap (figs. 25 and 26) presents nearly the same geom- 
etry to the flow as does the transverse gap at o? = 45° (fig. 24, y = 3.81 cm). However, 
the gap-corner flow reattachment observed for the longitudinal gap was not present on the 
transverse gap. These contrasting results are attributed to differences in the gap width 
and corner radii between the longitudinal and transverse gaps. As a result of the model 
fabrication process discussed previously, longitudinal gap width was approximately 
1.5 times as large as transverse gap width and longitudinal gap edge radius was larger 
than that of the transverse gap. The increased gap width and large exposed surface area 
at the longitudinal gap edge allow for diffusion of a shear layer into the gap and flow 
reattachment at the tile corner as opposed to the negligible separation and reattachment 
observed for the transverse gap. 

The variations of gap heat transfer with rotation angle for the simple transverse gap 
and the gap intersection are presented in figures 27 and 28, respectively. Array rotation 
angle has little effect upon the heating within the simple transverse gap (fig. 27), but tends 
to increase slightly the heating to the surface immediately downstream of the gap. For 
the gap intersection location (fig. 28) at the lower Reynolds numbers, array rotation again 
has minimal effect on heatir^ within the gap. However, at the highest Reynolds number 
(fig. 28(c)) array rotation decreases the gap heating level at the gap intersection. It 
should be noted that this heating decrease is a result of the previously discussed shift of 
the heating peak within the gap (figs. 22 to 24), and does not reflect a lessening of the 
heating load within the gap. Figure 28 also clearly demonstrates the large increase in 
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tile surface heating which occurs immediately downstream of a longitudinal-gap— 
transverse-gap intersection. All surface data of figure 28 are significantly higher than 
the corresponding data for the simple transverse gap (fig. 27). 

Pressure -Gradient Effects on Gap Heat Transfer 

Rotation of the curved array not only produced crossflow over the array, but also 
the desired pressure gradient within the gap. The magnitude of this pressure gradient 
was determined by differentiatii^ the curve fits of the gap pressure data of figure 21. In 
figure 29 measured heat transfer in the transverse gap at y = ±3.81 and ±11.43 cm is 
presented as a function of a nondimensional pressure-gradient parameter. This param- 
9P L 

eter — — is the local pressure gradient divided by the local pressure times some 

characteristic length. Using the full-scale tile dimension as the characteristic length 
(L = 15.24 cm), this parameter is physically the nondimenslonalized pressure drop for 
one tile length along the gap. No systematic effect of pressure gradient on gap heat 
transfer is evident. 


Shuttle Pressure -Gradient and Boundary -Layer Simulation 


In order to evaluate the full-scale shuttle pressure -gradient simulation obtained 
with the curved-array model, full-scale-vehicle spanwise pressxire distributions were 
obtained at two longitudinal stations. These pressure distributions were calculated, 
usii^ the method of reference 10, for a full-scale vehicle at Mach 10.0 and 30° angle of 
attack in an ideal gas with y = 1.12. The surface pressure distributions were numeri- 
cally differentiated to obtain the full-scale values of the pressure -gradient parameter 


^ Values of this parameter for the full-scale vehicle are presented in figure 30 

as a function of surface dimension S/Yj^g^x for the two longitudinal stations illustrated. 
The shaded area superimposed upon these data indicates the range of pressure -gradient 
parameter values obtained in the wind-tunnel tests. Simulation of full-scale-vehicle 
pressure-gradient levels is excellent with the exception of the wing leading-edge regions. 
Wing leading-edge TPS, however, will be a solid mater^l with no gaps present, therefore, 
pressure gradient is not of concern on that portion of the vehicle surface. 


Reference 5 indicates that although the boundary -layer edge Mach number and unit 
Reynolds numbers of the present tests are more extreme than expected in flight, the 
boundary -layer displacement and momentum thicknesses provide good simulation for a 
range of flight body-point— trajectory-point combinations. 
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CONCLUDING REMARKS 


An experimental investigation has been conducted to assess the effect of pressure 
gradient on the heat transfer to reusable surface insulation (RSI) tile-array gaps sub- 
merged in a thick, turbulent boundary layer. The experimental program consisted of 
heat-transfer and pressure measurements on a cxirved array of full-scale simulated 
RSI tiles in the tunnel-wall boundary layer of the Langley Research Center continuous 
flow hypersonic tunnel over a range of free-stream Reynolds numbers and flow angular- 
ities. The tile-array model was tested with gaps present, and with gaps eliminated to 
obtain smooth-surface reference data under conditions simulating full- scale -vehicle 
pressure gradients. In order to gain a thorough understandii^ of the boundary -layer flow 
in which these tests were conducted, pitot profile measurements of the tunnel-wall bound- 
ary layer and wall-pressure and heat-transfer measurements were made. 

Because of corner effects in the square cross section contoured nozzle, the wall 
boundary layer of the continuous flow hypersonic tunnel is highly three dimensional with 
significant transverse thickness variations. However, the measured flat-plate wall heat 
transfer correlated with momentum-thickness Reynolds number and boundary -layer 
displacement thickness. Heating to the curved tile array in this flow was more sensitive 
to transverse position than was the flat-plate heating. This sensitivity resulted from 
thinning of the boundary layer by the protrusion of the curved surface into the flow. A 
method was developed for prediction of smooth -curved -array heat transfer in this tunnel- 
wall boundary -layer flow. The method uses an estimated or measured surface pressure 
to perturb the measured flat-plate heating data, and accounts for "effective boundary - 
layer thinning due to the protrusion of a model into the boundary layer. Simulations of 
full-scale space shuttle transverse pressure gradient and boundary -layer thickness were 
good. 

The results of this investigation indicate the following: 

1. The level and distribution of pressure on the floor of the RSI tile -array gaps 
follow the trends of pressure on the external tile surface. 

2. Heat transfer to a forward-facing transverse-gap wall is significantly higher 
at the intersection of a longitudinal and transverse gap than for a simple longitudinal 
or transverse gap. 

3. Heat transfer to the tile surface immediately downstream of a transverse gap 

is higher than for the smooth surface with no gaps present. The increased surface heat- 
ing due to the gaps is particularly significant downstream of a lor^itudinal -gap- 
transverse -gap intersection. 


4. For a thick, turbulent boundary layer there is no systematic effect of transverse 
pressure gradient on tile-array-gap heat transfer. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 14, 1975. 
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TABLE L- THERMOCOUPLE AND PRESSURE ORIFICE LOCATIONS 


[Dimensions are in centimeteri0 


Thermocouple 


X 

y 

- 27.00 

0 

- 27.00 

2.54 

- 27.00 

5.08 

- 32.06 

0 

- 29.53 


- 24.47 


- 21,94 


- 11.82 

2.54 

- 11.81 

5.08 

- 11.79 

7.62 

- 11,77 

10.16 

- 11.75 

12.70 

- 16.85 

7.62 

- 14.32 

1 


Thermocouple 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 


1.30 3.04 
.73 3.04 
.16 3.04 
.16 --- 


Pressure 

orifice 
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(a) Model access position. 

Figure 2.- Model injection mechanism. 
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(b) Rotated to test position. 




Figure 3.- Schematic drawing of the RSI pressure-gradient model. 
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Figure 5.- RSI pressure -gradient model instrumentation locations. 
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Figure 6.- Flat-plate model mounted on injection strut. 
















Figure 12.- Transverse variation of boundary -layer velocity profile 

= 3.3 X 10® per meter. 




Figure 13.- Variation of displacement and momentum thicknesses with 

transverse wall position. 

















Figure 15.- Static pressure distribution on the flat plate. 































































□ Measured data 0 

Equation (13), with 

measured pressures 



y, cm 


Figure 20, - Comparison of measured and predicted heat transfer to the smoothed 
RSI tile-array model using measured pressure data. = 3.3 x 10^ per meter 





















meter. 




























a, degrees 

(a) = 1.6 X 10^ per meter. 

Figure 25.- Variation of heat transfer to the longitudinal gap with array 
rotation angle, x » -15.61 cm. 
















a, degrees 


(a) R«j = 1.6 X 1 q 6 per meter. 

Figure 26.- Variation of heat transfer to the longitudinal gap with array rotation angle. 

X = 8.02 cm. 





























